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Abstract 

We study properties of high energy factorisable gluon densities focusing on saturation 
effects. In particular we show that the property of saturation of unintegrated gluon 
density allows for introduction of thermodynamical entropy associated with production 
of gluons. This is due to the observation that saturation scale acts like a mass of a gluon 
which related with temperature via thermodynamical relations allows for calculations of 
entropy. We also show that obtained entropy behaves like multiplicity of produced gluons. 



1 Introduction 



Perturbative Quantum Chromodynamics (pQCD) at high energies generates scale called satu- 
ration scale Qs- The need for the existence of such a scale originates back to investigations of 
unitarity violation |T] by linear equations of pQCD [2l|3]. To resolve this problem, higher or- 
der perturbative corrections of nonlinear type within the Balitsky-Kovchegov, CGC/JIMWLK 
framework (i.e. nonlinear modifications to summation of logarithms of energy a"(lns)") were 
considered [H-[T^. These unitarity corrections have a clear physical meaning for they allow glu- 
ons to recombine what at some point balances gluons splitting. The saturation scale depends 
on energy and its existence prevents gluon densities from rapid growth. Existing data suggest 
that the phenomenon of saturation occurs in nature. The seminal example is provided by a 
discovery of the geometrical scaling in HERA data [15] and more recently by a presence of 
geometrical scaling in production of inclusive jets in the LHC data p^ITT] . 
Here, we would like to use the saturation property of gluons to derive an expression for amount 
of thermodynamical entropy associated with generation of saturation scale Qs which acts as 



effective mass of gluon during collision of hadrons. We show that such a construction is possible 
due to the relation of temperature to the saturation scale and the property of saturation scale 
which acts as an effective mass of gluon. Our main result are formulas (!25l) . (!26|l . for other 
approaches to entropy production we refer the reader to [T8H22]. We also show that the num- 
ber of gluons whose distribution is given by unintegrated gluon density (which in perturbative 
framework is obtained by summing via integral equations Feynman diagrams describing emis- 
sions of gluons) f l27|) is up to a constant factor equal to number of produced gluons therefore it 
can be directly linked to entropy of produced gluons. This is possible because of universality 
of saturation scale which is the only scale in the problem and which links those quantities in a 
unique way. 

The letter is organised as follows. In section 2 we obtain a formula for a gluon density in 
the adjoint representation of the color group and we show that saturation is reached at larger 
scales than in the fundamental representation. In section 3 we use the obtained gluon density 
to calculate the cross section for inclusive gluon production and we show that the cross section 
behaves like the saturation scale. In the 4'th section we derive formula for entropy associated 
with produced gluons and we relate it to number of gluons given by unintegrated gluon density. 
Finally we conclude that since saturation scale saturates itself the entropy of gluons coming 
from saturated region has to be bounded. 

2 Saturation scale in the adjoint color representation 

The high energy factorisable [23] unintegrated gluon density can be linked to a dipole amplitude 
for a color dipole scattering off a hadronic target via the following relation 

0(G),(Q)(x,/c) = J ci^bd^re-^^- V^iV(G),(Q)(r,6,x). (1) 

where the subscripts Q,G specify whether considered dipole is in the fundamental or in the 
adjoint color representation. The impact parameter at which a dipole collides with the hadron is 
b, the size of the dipole is r = |r| and k = |k| is its conjugated variable characterising transversal 
momentum of the gluon. The two dimensional vectors 6 = |b|, A; = |k|, r = |r| lay in the 
transversal plane to the collision axis. The variable x is related to the longitudinal momentum 
component of the gluon momentum. 

Due to the fact that in an electron-proton (e-p) collisions a dipole originates from a photon 
dissociation the associated gluon density is taken to be in the color fundamental representation 
while in hadron-hadron (A-A, p-A) collisions it originates most likely from gluon, therefore the 
associated gluon density with a gluon dipole is in the color adjoint representation. The relation 
between these two amplitudes in the large N^. limit is the following [2l] : 

NG{r,b,y) = 2NQ{r,b,y)-N'Q{r,b,y), (2) 

In our investigations of the gluon density in the adjoint representation, we use the one which 
can be obtained from the GBW model [27] for the dipole amplitude. This model although, 
quite simplistic, is a nice laboratory for studying physics of saturation. In particular, we find 



that the integrals needed in investigated in the next section inclusive gluon production can be 
performed exactly. The GBW amplitude reads: 



NQ{x,r,h) = d{h^-h) 



1 — exp 



(3) 



A/2 



where Qs{x) is the saturation scale and it is modelled to be Qs{x) =Qo (^)"'' where A, Xq, 
are free parameters and bo defines proton's radius. This amplitude saturates for large dipoles 
r and exhibits geometrical scaling which has been confirmed by data [T5],[27j. The 

unintegrated gluon density in the adjoint representation can be obtained in a straightforward 
way via the transformation ([1]). However, before we perform this transformation, let us recall 
the result for the gluon density in the fundamental representation: 



exp 



(4) 



where we assumed a cylindrical shape of the proton which integrated over the impact parameter 
b gives the transversal area A± = nR"^, where i? is a radius of a proton. Now motivated by the 
fact that this function exhibits a maximum as a function of k"^, see fig. ([T]) for fixed x, we define 
the saturation scale as momentum for which the (f)Q{x, k) has maximum as a function of k"^: 



Qs{x) = (9fc20Q(x, A;) = 



Applying this condition we obtain: 



dk2(f)Q{x,k) 



<pQ{x,kt){Q^,{x) 



k't) 



k^Qlix) 







(5) 



(6) 



This observation has been used in [28] to formulate the GBW model as a solution to a transport 
equation of a form: 

dY(l){Y,L) + XdL<f){Y,L,b) = (7) 

with initial condition 0(fc^,a;o) = fc^exp(— A;^) where Y = lnl/x and L = Ink"^. Now we calcu- 
late the saturation scale for gluons in the color adjoint representation. To see that the effective 
saturation scale for gluons in the adjoint representation is larger than in the fundamental rep- 
resentation we write the generic expression for the dipole cross section which saturates at large 
dipole size as: 



Using 



NQ{x,b,r) = 1 - 
in (|2]) and rearranging terms we obtain: 

NGix,b,r) = 1 _e-2/(-.nb) 



(9) 



where the factor 2 in the exponent is the origin of the difference between the saturation scale for 
gluon density in the adjoint representation and gluon density in the fundamental representation. 
To obtain the gluon density we use ([1]) together with ([2]) to obtain gluon density which reads: 



(l)Gix,k) 



exp 



2Qf{x) 



(10) 
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Figure 1: High energy factorisable gluon density. The maxima signalise existence of the saturation 
scales. Magenta continuous line represents the gluon density in the adjoint representation. The blue 
dashed line represents the gluon density coming from the dipole in the fundamental representation. 
Plot is made for x = 10~^. 

Applying our definition for saturation scale and calculating derivative we find: 



what means that the effective saturation scale for gluons in the adjoint representation is larger 
than in the fundamental representation and reads {Q1)g = 2((5^)q. 

3 Inclusive gluon production 

In this section we calculate the rapidity dependence of the cross section for the inclusive gluon 
production using the gluon density (fTOl) . Since we will use only the gluon density in the adjoint 
representation we skip the index G. Here, we consider a slightly asymmetric configuration near 
the central region, that is why we use the same parameterisation for the gluon densities but with 
two different saturation scales. In the following, we are interested in qualitative results, therefore 
for simphcity, we assume hadrons to be protons. Following |10], we use xi = (Qsi/-\/s)e^ and 
^2 = {Qs2/-\/s)e~^ where y/s is the energy of the collision. The formula for inclusive gluon 
production reads [DISlESlEniEn]: 



For recent theoretical studies of the inclusive gluon production with the NLO effects coming 
from the running coupling constant we refer the reader to [31]. For recent phenomenological 
works see [321439] . 

Using the expressions for the unintegrated gluon densities and performing the integral over 
(Pk, we obtain: 






(12) 
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Figure 2: Cross section for inclusive gluon production (normalised to larger of two saturation scales 
and square of the transversal area A± ^ S±) as a function of ratio of saturation scale. The blue line 
is the full cross section [14^ , the magenta line is the leading term, the green is the contribution from 
the sum of leading and subleasing term. 

Integrating over d'^pt and introducing minimal ptmin because the integral is logarithmically 
divergent, we obtain: 



2 



(14) 

Assuming, that the saturation scale in one of the colliding hadrons is the smallest scale in the 
problem Ptmin — Qsi<Qs2 we perform expansion in terms of Qsi/Qs2 in dll]) arriving at: 

da _ 2AICfQI, 2A\CfQ\^ 



dy Ti'^OLs Ti'^asQ'^ 



s2 



2 log + 5 + 27^) + O mjQl.f) (15) 



A comment regarding formula (12) is in order. In derivation of (12) it has been assumed that 
one of gluon densities is in a dilute regime, the more general formulation which would allow for 
taking into account interaction of two dense targets would most probably introduce regulator of 
the infrared divergence. However, such framework at present is not known and we regulate the 
divergence by the lowest physical scale in the problem. The result above shows in agreement 
with inilES] that da/dy behaves like the saturation scale making this observable particularly 
interesting for studying physics of dense partonic systems formation. On Fig. ([T]) we show 
results obtained from formulas ( !T^ and (fT5|) . We see that the full result is well approximated 
by the first term if the saturation scales are significantly different. 



4 Entropy 

Now we would like to relate the obtained cross section for the inclusive gluon production to 
produced entropy in the collision due to emergence of effective mass of gluons during the 
collision. We wish, however, to consider dense system where the following formula is not easily 
justified : 

/d'^7 d-^k 
-^[-f{z, k) In f{z, k) + .,(1 + sj{z, k)) ln(l + sj{z, k))] (16) 



where f{z,p) is phase space density of patrons/hadrons, Sj = —1 for fermions and Sj = 1 for 
bosons. We see the possibihty to address the question of entropy of dense system, following 
ideas of [111^2] where relation between saturation scale and temperature of thermally produced 
hadrons has been obtained within partonic picture of QCD and it reads: 

T = (17) 

where x stands for a longitudinal momentum fraction of hadron carried by a parton. The for- 
mula above is a direct consequence of the Unruh effect [43j within QCD. The Unruh effect [13] 
states that accelerated observer in its rest frame feels thermal radiation of temperature T = 
In the construction leading to QCD application, a decelerating off shell gluon with distribu- 
tion in transversal momentum given by unintegrated gluon distribution is considered. Due to 
deceleration, a horizon is formed tunnelling through which produces a thermal ensemble with 
a density matrix of the Maxwell-Boltzmann type. To determine the deceleration the authors 
of IH] use the Wong equations ^Ij for color particle moving in a homogeneous chromoelectric 
field. Solutions of this equations give a result that \a\ = Qs which leads to eq. ( |T7I) . Fur- 
thermore, it is well known that saturation scale effectively acts as a regulator of the infrared 
divergent behaviour of the unintegrated gluon density. At the lowest order in QCD (in phase 
space region where the total energy is larger than any other scale) it is: 

asC*F 1 



This quantity has a structure of propagator and we see that at small values of k it becomes di- 
vergent. Higher order perturbative corrections coming from BFKL resumation do not solve this 
problem. This low k behaviour of BFKL equation is interpreted as leading to infinitely many 
gluons with small k and also to power like growth of parton density with energy what violates 
unitarity bounds. Saturation due to nonlinear effects modifies the behaviour of distribution of 
gluons and in particular, since the unintegrated gluon density has a maximum, it sets the most 
probable momentum to be of the order of Qs, what can be read off from formula (ITT]) . This we 
interpret as giving by nonlinear effects (here we study saturation model but saturation is also 
studied within nonlinear QCD evolution equations [SIE]) emergent mass (gauge invariance still 
holds) to gluon which is proportional to the saturation scale and which renders the propagator 
to be finite. 

We are also motivated by studies of QCD by means of nonperturabative methods where it is 
argued that gluon mass consistent with gauge invariance might be dynamically generated due 
to nonlinearities (see [15] and references therein) . The regulated divergency of propagators due 
to generation of effective mass Qs{x) of virtual gluons has also impact on real gluons since they 
acquire mass. We can therefore write for single real gluon : 

Mcix) = Qsix) (19) 



This observations allows us to define entropy of gluons, because the existence of temperature 
and existence of effective mass of gluons allows for introduction of the thermodynamical entropy 



(coarse grained entropy). Using the laws of thermodynamics and assuming small change of the 
volume associated with system of gluons we obtain: 



dE = TdS (20) 

and setting dE = dM gives: 

dM = TdS (21) 

where the mass is understood to be the emergent mass of system of gluons carrying small 
longitudinal momentum fraction of the hadron and is related to mass of single gluon via: 

M{x) = Ng{x)Mg{x) (22) 
where Ng{x) is (dimensionless) number of gluons defined as: 

,r / N d.N 1 da 

Ng{x) = — = — — (23) 
dy o± dy 

where dN/dy is multiplicity of produced gluons. Now we use eq (fT7|) what allows us to link the 
saturation scale to produced entropy: 

d[NG{x)MG{x)] = ^^dS (24) 



Thus we obtain: 



S=^-^Ql{x) + So (25) 



and equivalently 

S = 37T [Ng{x)+ Ngo] (26) 

where 5*0 = ^'^^'^ /^^ is a state of the lowest entropy or the minimal number of gluons. 
Obtained entropy can be also expressed in terms of number of gluons whose distribution is 
given by unintegrated gluon density. The number of gluons which build up unintegrated gluon 
density is: 

ug^x) = - I d^rd^k^{x,k,r) = - [ d^k (f){x, k'^) (27) 



TT J IT _ 

where d^r = d^bdl and r = (/,b), / is longitudinal dimension of the proton, b is already 
introduced impact parameter, ^{x,k,r) is understood to be gluon density unintegrated in r, 
k, X. For simplicity we assume that longitudinal space distribution of gluons factorises and 
is given by 6{l). This assumption is justified since the nontrivial dynamics takes place in the 
transversal space. Therefore we obtain: 

ncix) = ^Q^(x) (28) 

Now we can expressed entropy of produced gluons by number of gluons given by unintegrated 
gluon density: 

^ = 127rnG(x) + 37rA/'Go (29) 



This can be used perhaps to define via fl2^ entropy of gluons described by unintegrated gluon 
density, however we postpone this task for future studies. Although, we used GBW model to 
arrive at formula (!25|) our result is model independent in saturated region since GBW model is 
consistent with results of I5l[6l|26lll7lll8] where gluon density in saturated region behaves like: 



The results we obtained above and in particular (12^ are valid when fi where fi is an 

unknown constant, we choose this constant to be Aqcd- In order to obtain a formulation valid 
also for small temperatures one has to go along lines of the Nernst theorem, going along these 
lines, one can expect that the lowest possible entropy will be zero what will correspond to 
situation where coherence of a proton [IB] is not affected by interactions. Furthermore, we see 
in agreement with intuition that the entropy is proportional to number of gluons and depends 
linearly on the hadron size i.e. doubling the hadron size the entropy is doubled. 
Here we would like to stress that our construction of entropy given by (1261) is quite general. The 
crucial element is the existence of saturation which links the temperature and the mass of system 
of gluons allowing for introduction of thermodynamical relations between these quantities. 
There is one more consequence of the dependence of entropy on the saturation scale. It gives 
upper bound on the amount of produced entropy associated with gluons confined in saturated 
region of phase space what follows from our result combined with result of |l9l|50]. In these 
investigations authors found that at extreme energies saturation scale eventually saturates what 
in combination with ( 12^ gives: 



Entropy however can still be produced outside of saturated region even if saturation scale will 
saturate itself. This we can guess from approaches which take care for more realistic behaviour 
of gluon density outside of saturation region. Namely outside of the saturation region i.e. 
k > Qs the gluon density from GBW model falls too fast (due to exponential behaviour it can 
be integrated to infinity) while BK or McLerran-Venugopalan |17] model gives hard perturbative 
behaviour of type Ql/k"^ which gives logarithmic divergent contribution after integrating what 
will give contribution to entropy unbounded by saturation of the saturation scale. 

5 Conclusions 

In this letter we have studied properties of the gluon distribution possessing the feature of 
saturation. We defined the saturation scale mathematically as a maximum of the gluon density 
f lTT]) . We have calculated the cross section for the inclusive gluon production and have shown 
that it might be approximated by saturation scale characterising gluon density of less saturated 
hadron. Furthermore, we introduced saturation scale related entropy of produced gluons and 
expressed it in terms of multiplicity of produced gluons. We have shown that this entropy be- 
haves like number of gluons whose distribution is given by unintegrated gluon density. We also 
concluded that since saturation scale ultimately saturates itself, there should exist an upper 




(30) 




(31) 



bound on entropy of gluons coming from saturated part of gluon distribution. The result we 
obtained is of course a contribution to the entropy coming just from emergence of saturation 
scale and effective mass of gluons in perturbative regime. 
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